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Abstract 

With a well-motivated extension of higher order holonomy corrections, the quantum theory of 
loop quantum cosmology (LQC) for the k = Friedmann-Robertson- Walker model (with a free 
massless scalar) is rigorously formulated. The analytical investigation reveals that, regardless of 
the order of holonomy corrections and for any arbitrary states, the matter density remains finite, 
bounded from above by an upper bound, which equals the critical density obtained at the level of 
heuristic effective dynamics. Particularly, with all orders of corrections included, the dynamical 
evolution is shown to follow the bouncing scenario in which two Wheeler-DeWitt (WDW) solutions 
(expanding and contracting) are bridged together through the quantum bounce. These observations 
provide further evidence that the quantum bounce is essentially a consequence of the intrinsic 
discreteness of LQC and LQC is fundamentally different from the WDW theory. Meanwhile, the 
possibility is also explored that the higher order holonomy corrections can be interpreted as a result 
of admitting generic SU(2) representations for the Hamiltonian constraint operators. 
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I. INTRODUCTION 



In previous years, the status of loop quantum cosmology (LQC) has undergone lively 
progress and become an active area of research. Specifically, with the inclusion of a free 
massless scalar field, the comprehensive formulation of LQC in the k = Friedmann- 
Robert son- Walker (FRW) (i.e. spatially flat and isotropic) model has been constructed 
in detail, providing a solid foundation for the quantum theory and showing that the big 
bang singularity is resolved and replaced by the quantum bounce, which bridges the present 
universe with a preexisting one [l, Q, S]- Similar results are also affirmed for a variety of 
extended models 

However, it remains questionable whether the quantum bounce indeed results intimately 
from the quantum nature of Riemannian geometry of loop quantum gravity (LQG), as the 
same result can be readily obtained even at the level of heuristic effective dynamics without 
appealing to the sophisticated features of LQC 



In response to this criticism, a simplified but exactly soluble model of LQC has been 
studied and used to prove that resolution of the classical singularity and occurrence of the 
quantum bounce are robust [ll[ . Furthermore, the study of [ll[ brings out the precise sense 
in which the Wheeler-DeWitt (WDW) theory approximates LQC and the sense in which this 
approximation fails, thereby showing that LQC is intrinsically discrete and the underlying 
discreteness is essential for the quantum bounce. 

To add further evidence for the loopy nature (i.e. intrinsic discreteness) of the quantum 
bounce, a new avenue has been suggested in jl2j to consider a well-motivated extension of 
higher order holonomy corrections. At the level of heuristic effective dynamics, the investi- 
gation of 12J reveals that the key features of the bouncing scenario are generic, regardless of 
the inclusion of higher order corrections. This strongly urges one to formulate the quantum 
theory of LQC with the higher order holonomy corrections and attest the observations of 



[12| on a firmer ground. In this paper, by applying the techniques introduced in [TJJ to 
construct the physical Hilbert space and identify a complete family of Dirac observables, the 
quantum theory of LQC is rigorously formulated for any order n of holonomy corrections 
(including the case of n = oo). The detailed investigation shows that, for any arbitrary 
physical states (not only restricted to the states which are semiclassical at late times) and 
at any order n of holonomy corrections (including n = oo), the quantum bounce is generic 
and the expectation value of the matter density is bounded from above by an absolute upper 
bound, which equals the critical density obtained in the heuristic analysis of [l2j. By the 
same notion of "singularity resolution" suggested in [ill], the classical singularity is said to 
be resolved. With the extension of higher order holonomy corrections, we affirm to a broader 
extent the robustness of key features of LQC previously shown in [111], the model of which 
is now viewed as a special case at the lowest order (n = 0) of holonomy corrections. 

Meanwhile, the higher order holonomy corrections might be understood as a result of 
using generic SU(2) representations (labeled by spin j, not restricted to j = 1/2) for the 
holonomies in the Hamiltonian constraint operator. To explore this possibility, we construct 
the Hamiltonian constraint operator in generic j representations in a fashion slightly modi- 
fied from that of [13|] and then investigate the linear sum of the Hamiltonian operators over 
generic j. The investigation is not conclusive but encouraging and reminiscent of the issues 
regarding large j in some spin-foam models, calling for further research on the j ambiguity 
from a new perspective. 

This paper is organized as follows. After a brief overview on the motivations in Sec. [Til 
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the main text is divided into two major parts. In Sec. II I II as the first part, we study the 
Hamiltonian constraint operator in generic j representations and their linear superposition. 
In Sec. |IV]as the second part, the quantum theory of LQC with higher order holonomy cor- 
rections is formulated and used to investigate the physical dynamics. Finally, the conclusions 
are summarized and discussed in Sec. |V] 



II. MOTIVATIONS 

In the presence of a free massless scalar field 4>{x,t) = <p(t), the classical Hamiltonian 
constraint for the k = FRW model is given by 

C = C gr „ v + a,,.„ = -^eW+^ (2.1) 

in terms of the Ashtekar variables c and p, the conjugate momentum p^ of <p, the lapse 
function N, and the Barbero-Immirzi parameter 7. At the heuristic level, to impose the 
loop quantum corrections of LQC, we simply take the prescription of "holonomization" to 
replace c with 

c — ^ (2.2) 
H 

by introducing a discreteness variable Ji. The heuristic effective dynamics is then solved as 
if the dynamics was classical but governed by the new "holonomized" Hamiltonian, which 
reads as 

3iV sin 2 /2c r-- Np 2 



C »^8^^bl + ^. (2.3) 

The bouncing scenario can be easily obtained at the level of heuristic effective dynamics 
without invoking the sophisticated features of LQC jl, 0, [Io| , therefore calling into question 
whether the quantum bounce results from the very quantum nature of LQC In particular, 
with the "improved" scheme (14. lj) imposed for jl, the modified Hamiltonian constraint = 
immediately sets an upper bound for the matter density: 

n 2 R 

sin fic < 3ppi, (2.4) 



H *' 2M 3 8ttG 7 2 A 
where the Planckian density is defined as 

Ppi := (8ttG 7 2 A)- 1 . (2.5) 

Apparently, without going into the detailed construction of LQC at all, it is anticipated 
that the matter density is bounded from above and thus the quantum bounce occurs. One 
might then argue that the boundedness of has little to do with the fundamental structure 
of LQC but merely an ad hoc consequence of the formal modification of (12. 2p . In this sense, 
the quantum bounce seems not really "loopy" enough. 

In regard to this issue, a well-motivated extension of higher order holonomy corrections 



has been suggested in [12| and used to test whether the quantum bounce is essentially a 
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consequence of the loopy nature (i.e. intrinsic discreteness) of LQC. By evoking the Taylor 
series: 

(2*0 ! 
2 2k (k\) 2 (2k + l] 



fc=0 

for — 1 < x < 1, the nth order holonomized connection variable is defined as 

c (n) - = I V" 2^ fsin uc) 2fc+1 f2 7) 

c ^ ' u^2 2fc (£;!) 2 (2fc + l) lS m/iCj ' 10 

which can be made arbitrarily close to c (as n — > oo) but remains a function of the holonomy 
sin /2c and the discreteness variable /2. To implement the underlying structure of LQC by 
replacing c with holonomies, can be used as an improved version of (12.21) . which now 

reads as c^~ Q \ In place of (12.31) . the Hamiltonian constraint with holonomy corrections up 
to the nth order is then designated as 

C<"> = Cg v + C„„ := - (c<"») W+ ^s, (2-8) 

which supposedly reflects the intrinsic discreteness of LQC in a more elaborate fashion. 

One might suspect that if we include corrections of all orders (n = oo), the quantum 
theory of C {n=oo) will lead to the same result as the WDW theory, as formally we have 
c l°°' ) = c - This should not be the case, because the elementary variables are still sin /2c and 
p, instead of c and p, even in the limit n — > oo and, therefore, the striking difference between 



LQC and the WDW theory as emphasized in 11] should persist 



The analysis of [12] at the level of heuristic effective dynamics indeed shows that the 
bouncing scenario is generic, irrespective of the order of corrections, and the matter density 
remains finite, bounded from above by an upper bound in the regime of the Planckian 
density, even if all orders of corrections are included. In addition to [ll|, this observation 
adds further evidence that the intrinsic discreteness of LQC inherently accounts for the 
quantum bounce and thus LQC is fundamentally different from the WDW theory. It also 
encourages one to construct the quantum theory of LQC with the higher order holonomy 
corrections. 

In Sec. \1V\ we formulate the quantum theory of LQC associated with the nth order 
Hamiltonian constraint: 

fil{n) _ r i{n) , r i ._ 3 ( (n)\2 2 , P<t> Icy q \ 

° - °grav + °matt ■— g^Q^2 \ C h ) P + y ' \ l - J ) 

which is rescaled from (12.81) with N = \p^ 2 ■ Not only can LQC based on (12. 9p be rigorously 
constructed but it also confirms the bouncing scenario to be generic, regardless of the order 
n, with the matter density bounded from above by an absolute upper bound, which equals 
the critical density obtained at the heuristic level in [l3|. 

The higher order holonomy corrections corresponding to higher powers of sin ]xc are most 
likely to be interpreted as the imprint of generic j representations for holonomies in the 
Hamiltonian constraint operator. Sec. II I II is devoted to investigate this possibility. The 
aim is to show that, with proper regularization, the linear superposition of the Hamiltonian 
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operators in generic j representations might conspire to give the operator of C^ av °°^ or 
C^av , i.e. given with appropriate coefficients Cj, we might have 



E 



3 ^grav — '-'grav \^- Lyj ) 



or 



E 



= (7'(°°) Co 11) 



where ^C' gra:v is the gravitational part of the (rescaled) Hamiltonian operator in the generic 

j representation and (%rav is the operator corresponding to C'J^v 00 ^ ■ This evokes the idea in 
|l4 | that the extension of the Hamiltonian constraint operator to generic j representations 
is motivated by the spacetime covariant formulation (such as spin-foam theory) of quantum 
gravity. In particular, the requirement of crossing symmetry as manifested in the covariant 
formulation naturally leads one to consider the Hamiltonian operators in j, which add links 
of arbitrary colors of j. A linear combination of such operators could actually exist and 
define a crossing-symmetric quantization of the Hamiltonian constraint of general relativity. 



On the other hand, the general considerations in [13| for LQC and in [15[ for LQG both 



disfavor the Hamiltonian operator in gen eric j representations other than j = 1/2. The 



fact that the problem suggested in [13l . Il5l | could be avoided is commented on at the end of 
Sec. UVDl 

Remark: 

There are many other avenues of adding higher order holonomy corrections in the liter- 
ature which are motivated differently and yield distinct dynamical behavior. In particular, 
the difference between the consideration here and that of [l6|, TtJ is commented on in (l2| . 



III. HAMILTONIAN OPERATORS WITH HIGHER ORDER HOLONOMY COR- 
RECTIONS 



This section presents the first half part of the paper. In Sec. IIII A\ we briefly review the 
notation and background knowledge for the classical theory (see [2J for more details). In 



Sec. IIIIBt following the treatment in 13|, we construct the (rescaled) Hamiltonian operator 
of LQC in a generic j representation with some simplifications. Finally, Sec. IIII Cl investigates 
the linear sum of the Hamiltonian operators over generic j. Our goal is not to reach a 
definitive conclusion that a specific linear superposition of the Hamiltonian operators can 
give rise to the operator associated with C'^S/ as suggested in f)2.1ip . but instead we try to 
justify this possibility and hopefully to inspire further research on the issues of j ambiguity 
from this point of view. 



A. Classical Hamiltonian constraint 

In the standard Hamiltonian treatment of the cosmological model, the spacetime M is 
foliated as M = S x R with E being the homogeneous spacelike slice. In the k = FRW 
model, because E is noncompact, one has to introduce a finite-sized elementary cell V on 
E and restricts all spatial integrations to V. In the comoving coordinates x = (x, y, z) in 
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which the spacetime metric reads as ds 2 = g flu dx fJ- dx u = —dt 2 + a(t) 2 dx 2 , V is chosen to be a 
cubic cell with fixed coordinate length L on each edge. With the symmetries of homogeneity 
and isotropy imposed, the gravitational phase space variables (the connections A a l and the 
density weighted triads E can be expressed as 

Aj = cVo'^ul, (3.1a) 
E\ = p^V ' lo el (3.1b) 



Here, °e" and %j l a are a set of fiducial triads and cotriads {°e°; is the inverse of °uj l a ) that are 
left invariant by the action of the Killing fields of £ and adapted to be aligned with the 
edges of V (i.e. oc Sf and 1o l a oc 5 l a ). The fiducial 3- metric is given by °q a b = S^ujIIjI 
and V D is the fiducial volume defined as the volume of V measured with respect to °q a b, 
i.e. V := j v d?Xy/% = j$ j$ j$ d?Xy/°q = L?^/°q with °q := det°g a b. Consequently, the 
symmetry-reduced phase space is described by the pair of c and p, which satisfy the canonical 
relation: 

{c,p} = ^— ■ (3.2) 

The absolute value of p gives the physical area of each rectangular surface of V while the 
sign of p corresponds to the spatial orientation. 1 

In terms of c and p, the gravitational part of the Hamiltonian constraint can be written 

as 

~a ~b 3N 
k 



C - = "IfeW S/ xNe ^' kF ^' g "' = -ftS^^W' (3-3) 
where 

e:= |detE| 1/2 = L" 3 |j9| 3/2 (3.4) 

and N is the lapse function. Additionally, with the inclusion of a free massless scalar field, 
the total Hamiltonian constraint also contains the matter part given by 

CUt = / d 3 xNe^- = (3.5) 
Jv 2 2 \p\ s/2 

where p^ := |p| 3,/2 is the conjugate momentum of the homogeneous field 4>(x, t) = <p{t) with 

{4>,P*} = I- (3-6) 

The classical dynamics is governed by the total Hamiltonian constraint C = C grav + C matt 
as given in 02. ip . 



B. Loop quantization in generic j representations 

In the passage to quantum theory, the quantization procedure of LQC is more intricate 
than that of the WDW theory. This is because the kinematical Hilbert space of LQC is 



On the other hand, the physical interpretation of c depends on the dynamics. For classical dynamics, we 
have c = jLa(t), i.e. the change rate of the physical length of the edges of V (times 7) with respect to 
the proper time t. 
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in the "polymer representation" for p, instead of the standard Schrodinger representation. 
Consequently, the operator corresponding to c does not exist and additionally it is not 
densely defined if the inverse volume function | 3 / 2 - ls na i' V ely quantized as the operator 
with eigenvalues equal to the inverse of the volume eigenvalues. Thus, to construct the 
Hamiltonian constraint operator, we have to express the classical constraint in terms of the 
triad variable p and the holonomy h k , both of which have direct quantum analogs. 

Following the close similarity of the regularization procedure a la Thiemann's trick used 
in the full theory [l8j], we first recast the term involving the inverse triad e _1 as 

e-W k E a 3 E\ = ]T ^J^- <*» °e] °e\ Tr [hf\{hf)\ V}r) , (3.7) 

ffiL 



where 



h\ p) := V exp / r i A a i dx a = exp(j2cn) (3.8) 
Jo 

is the holonomy along the edge of coordinate length JxL aligned with the direction of °efd a , 
V = \p\ 3 ^ 2 is the physical volume of V, and Tj are the SU{2) generators with [n, Tj] = ^ij k T k 
(2iTi = o~i are the Pauli matrices in the standard convention). Note that this identity holds 
for any choice of p,, even when it is allowed to be a function of p. 

Next, invoking the standard techniques in gauge theories, Stokes's theorem allows us to 
express the field strength component F* b in terms of holonomies as 



Fab « -2Tr 



hff 



'q~ 1/3 ■ t 

^ O, ,7 Q ,k 



J Jx 2 L 2 



KH, (3-9) 



where 

h { &] k -.= hfhf{^)-\hf)- 1 (3.10) 

is the holonomy around the boundary of the square which is perpendicular to the di- 
rection of °efd a and of coordinate length fiL on each edge. In the context of homogeneous 
models, we have F l ah = d[ a A b f + e l j k Aj A h k = e l j k Aj A h k \ thus, Stokes's theorem is not es- 
sential in (13.91) and we can simply express TjAj °e" = °q~ l l % lim^ {hf^ — {hf^)~ l ~) / (4/iL), 
which leads to an alternative expression: 

F: b « -2Tr ([h?» - (hfY\h™ - (ZifW) J^M^t (3- 11 ) 

where the extra factor 2 in h[ 2 ^ is adopted for later convenience as will be seen. In the 
continuous limit — > 0, both (I3.9P and (13. lip become exact; however, in the theory of LQC, 
this limit does not exist and we should keep Jx finite as the "regulator" which reflects the 
"fundamental discreteness" of the full theory of LQG as an imprint on the reduced theory. 
To quantize C grav as an operator, keeping Jx finite, we depart from ( 13. 7ft with either ( 13.9ft or 
( 13. lip . Moreover, when the sum of higher j representations is considered in the quantization 
procedure, it turns out (13. lip is more controllable than (13.91) and thus is preferred. 2 



2 



The Hamiltonian operator based on ()3.1ip instead of (j3.9p has been investigated in [l9( in j = 1/2 repre- 
sentation. It shows that the alternative quantization gives rise to the operator C grav = 4 sin ^ A sin ^p, 
the dynamics of which is only slightly different from that of the operator C gra v = sin /2c A sin /2c obtained 
in the standard approach (j3.9p with A being the operator quantized from (j3.7p . In fact, if an extra factor 
2 in h[ 2 ^ is also adopted as in (j3.1ip . the former is identical to the latter. 
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To obtain the quantum theory, the holonomies h\ are promoted to ^fif in the j 
representation, V to V, and the Poisson bracket to a commutator. In the j representation of 
the SU(2) group, the Lie algebra generators Tj are represented as (2j + 1) x (2j + 1) matrices 
v> which satisfy 

Tr {Wrphj) = + l)(2j + 1) (3.12) 

By ( 13. 7p and (13.91) . the classical Hamiltonian (I3.3P with N — 1 then leads to the Hamiltonian 
operator in the generic j representation: 

3z 



x ^>*Tr (^hf Xi) hf ^hf^^hf ~^hf [^hf ~\v]) . (3.13) 



The Hamiltonian operator in generic j representations has been studied in depth in [13 
The matrix elements of holonomies in the j representation are given by 

2j-\m+n\ , , v , _ 2 ._ 



fU)hf\ = Tmn V" i-^co^/2 3 S sh^w/2, (3.14a) 

\ / mn ' * I mns 

s=\m—n\,\m—n\+2,... 

i-Y-rn+s . . 2 j-s s 

(^Kf) = T mn V ^— ^ cos/ic/2 sin/ic/2, (3.14b) 



s=|m— n| ,|m— n|+2,' 



OUf ) = e ^5 mn , (3.14c) 



where m, n G {— j, — j + 1, . . . , j — 1, j}, the index s in the sum increases by an increment of 
2, and T mn and Y mns are constants given by 

T mn = j{j + m )\{j-m)\{j + n)\{j-n)\, (3.15a) 
K TOns = ( j + ~(m + n - s) ) ! ( j - -{m + n + s) 



x ( -(m-n + s) 1! ( -(n-m+s) )!. (3.15b) 



The resulting Hamiltonian operator reads as 

Cf)c = 

grav " (87rG)^7 3 j(j + l)(2j + l)/i 3 



x V] (j) Z^'cos fxc/2 sin/ic/2 e i77l ^ c Ve-^ c , (3.16) 

m=—j s'=l,2,— 



where the constants ^'Z^ can be expressed in terms of T mn and Y mns (see 13j). 

Similarly, for the matter part of the Hamiltonian constraint, to deal with the inverse 
volume V~ l = |p|~ 3//2 in (I3.5p . Thiemann's trick is used again to define the inverse volume 
operator as 



(j)y-i 



3i 



87rGhjj(j + l)(2j + l)fl 



3 

2(l-£) 



(3.17) 



S 



Here arise two ambiguities, labeled by a half-integer j for different choices of i representations 
and a real number I (0 < I < 1) [20(. The general considerations in 13j, [l5| urge one to set 
j = 1/2. 3 For £, a general criterion is not available and I = 1/2 and i = 3/4 are most used 
in the literature. 

Unfortunately, the formula (13.161) does not yield the desirable expression which would 
readily match (12.111) when the operators ^^C gr3iV of generic j are properly summed. In 
addition, the inverse volume operator in (13.171) gives rise to further complications. To 
achieve our goal of having (I2.10p . we make two simplifications: First, we neglect the problems 
involving the inverse triad by choosing an appropriate lapse function N before quantization 
and thus obviate Thiemann's trick; second, instead of (13. 9p . the alternative quantization 
scheme ( 13. lip is adopted for the field strength component. 

In the full theory of LQG, it is not possible to get rid of the involvement of the inverse 
triad by gauge fixing N to eliminate e -1 , because e is a weight-one scalar density while N is 
a scalar with respect to diffeomorphisms on S. However, in the context of the k = FRW 
model, thanks to homogeneity, we have fl3.4[) and thus the difficulty due to the inverse triad 
can be readily avoided by choosing N = |p| 3,/2 . 4 With this choice of N, in parallel to (13. 3p 
and (13.51) , the rescaled Hamiltonian C = Cg rav + C' mautt reads as 

and 



CU = / (/3 ^i = Y' ( 3 - 19 ) 



3 



V 



as given in (12. 9p . With this gauge fixing imposed before quantization, the quantization 
strategy is much easier and does not need Thiemann's trick anymore. 

Furthermore, the second simplification is to adopt (13. lip instead of (13.91) . In parallel to 
( I3.13p . taking (13.111) into (I3.18P yields the (rescaled) Hamiltonian operator in the generic j 
representation given as 

ti)(j> 



srav 87rG 7 2 j(j + l)(2j + l)(4/i) 2 

x ^ fc Tr ([ Whf> - (^hfY 1 , WfcJW - (WfcJW)- 1 ] "^) \p\ 2 (3.20a) 

ijk 



^G 7 2 j(j + l)(2j + l)(4/i) 2 



xTr 



^h? ] ~ ^ti^r 1 ) [ U W ] ~ (^WY 1 ) U) r 3 \P\ Z , (3.20b) 



Since we will consider the linear sum of the Hamiltonian operators over all values of j, it is not clear 
whether one should still stick with j = 1/2 for the inver volume operator. This problem will not bother 
us, as we will rescale the Hamiltonian with an appropriate N before quantization to avoid the Thiemann's 
trick. See the remark in the end of Sec. IIII CI for more comments. 
4 Note that the choice of N = \p\ 3 ^ 2 gives a scalar since \p\ is the physical area of the faces of V and 
independent of the coordinates; on the other hand, e = L~ 3 \p\ 3 ^ 2 yields a weight-one scalar density 
because of the coordinate-dependent factor L~ 3 . 
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where the fact that the operator is SU(2) invariant has been used to arrive at (13.20bj) . Using 

( BHD , we have 



Tr 



r 3 



j 2j—\m+n\ 
m,n=—j si,S2=\m—n\,\m—n\+2,- 



T-m.rU-'rvrn. ( ^ 



Y Y 

1 mns\ 1 nms2 



x[l-(-l) Sl ][l-(-l) S2 ]oo~s^ 4i Sl S2 shT^ S1+S2 



=: ^'W s cos /ic J sin /ic 



2j 



s=2,4,- 



fc=l,2,- 



01 



-2A; 



Uk sin /ic 



(3.21a) 
(3.21b) 



where the index s is simply s = si + S2 and we have used the fact that the summand in the 
sum is nonvanishing only if s\ and S2 are both odd to arrive at (13.21bj) . This leads to 



00 A 



18 



2j 



grav 



8ttG'7 2 j(J + 1)(2j + 1)(4 / u) 2 



(A 



wy afcSin/ic 



2/v 



(3.22) 



vfc=l,2,- 



Explicitly, we list (j - ) C'g rav for the leading values of j: 



(1/2) C" 



grav 



3 



r^sin/xc 



grav 



sin //c — sin /ic ) \p\ 



(3/2) 



grav 



3 / 2 12 ; 6 i . ,, 

Ism/ic — sin/ic +-sm/xc 1 |p| 



87rG'7 2 y u 2 V 



(3.23a) 
(3.23b) 
(3.23c) 



(2) c: 



grav 



3 /- — r~2 21- — ^4 24- — ^6 8 r~a\ . l2 , 

■ „ 9 _ 9 sin/xc — — sin/xc + — sin [ic — -smite ) \p\ .(3.23d) 
87rG7 z /i z \ 5 5 5 



Note that the leading coefficient is given by ^hk=i = 8j(j + 1) (2j + l)/3, which is in accord 
with the fact that, for any generic j, ^C' gra:v returns to the classical counterpart in the formal 
limit — > 0. 

Remark: 

It has been hotly debated whether it is legitimate to rescale the Hamiltonian by the lapse 
function N before quantization. This question is in fact asked in two respects: (i) for the 
theory of LQC derived from LQG; and (ii) within the confines of LQC While the quest for 
a systematic formulation to derive LQC from LQG remains wide open, it is logically viable 
to obviate the issues of the inverse triad by gauge fixing N within the confines of LQC. 
One might still argue that, without the corrections by Thiemann's trick, the resulting LQC 
would no longer have the appealing feature of singularity resolution, as it has been shown 
that the state in the kinematical Hilbert space associated with the classical singularity is 
decoupled in the evolution equation as a result of the loop corrections on the inverse triad 
(2l| . However, decoupling of the singular state is neither a necessary nor a sufficient condition 
for singularity resolution. To provide a satisfactory notion of "singularity resolution", one 
should have available the physical Hilbert space and a complete family of Dirac observables. 
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This issue is discussed in depth in [11], and therein an applicable and satisfactory notion 
is suggested and used to show that the singularity is resolved and the corrections on the 
inverse triad are not essential. We will adopt the same notion of singularity resolution and 
reach the same conclusion for the quantum theory with higher order holonomy corrections 
in Sec. HVTH 

C. Linear superposition of the Hamiltonian operators in generic j 

As shown in (13.221) . ^C' gT!W is given as the sum of even powers of sin Jic up to sin 4j Jic. As 
the right-hand side of (12.111) also consists of even powers of sin Jic, the coefficients Cj can be 
obtained through the truncated version of (12.111) as follows. For a given n, we can uniquely 
determine the set of constants Cj such that 

n/2 n 

£ ^^a^-^^E^^"^ 2 -^, (3-24) 

3=1/2,1,- 1 ^ k=l 

where bk are the coefficients of the Taylor series: 

oo 

(sin- 1 ^) 2 = ^b k x 2k (3.25) 
fc=i 

— 1 < x < 1 and the operator Cgrll is defined as the truncated version of C&S) (up to a 
factor ordering). That is, for any arbitrary n, we can always find a linear superposition 
of ^C'grav to match up the Hamiltonian operator C'g^J with the terms of higher than 2n 
powers of sin Jic all removed. 5 

Equation (I3.2ip gives a convenient algorithm to compute the coefficients ^W, which 
are then substituted to (I3.24p to compute c^. The computation by computer shows that, 

undesirably, Cj do not converge to a set of constants but rather diverge rapidly in magnitude 
as n increases. Interestingly, a closer examination on the computed data suggests that the 
growth of Cj 1 ^ yields a peculiar profile of asymptotic behavior: 

cf ] » Cjne n/2 , Cj : constants (3.26) 



5 Alternatively, instead of C gr av , we can use Cgrav as a truncated form of Cgr^/ and replace ()3.24j) with 

n+l/2 

E (n) (j)A/ _ nl{n) 
j w grav ^grav ■ 

J=l/2,1,- 

However, C* g rav, the operator of Cgial as defined in (|2.9[) . does not precisely strip off C'^rS) the terms of 
higher than An + 2 powers of sin/Ic. At the limit n — > oo, the alternative treatment nevertheless yields 
the similar profile of asymptotic behavior for Cj] the computer computation in this approach indeed 
suggests Cj ~ 2cj n e n+1 / 2 in parallel to (|3.26[) for very large n, but the computation is far less efficient 
as the convergence is much slower. 
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when n is large enough, meaning that the relative weights of cj 1 nevertheless tend to be 

constant. As shown in Fig. (TJ the coefficients Cj rescaled by (ne n / 2 ) -1 tend to converge 
to a set of constants Cj and the relative weights Cj are appreciable only for moderate values 
°f 3 (j ^ 5) but diminish rapidly for large j. 6 This asymptotic behavior suggests that, 
comparatively, only the contributions of moderate j are manifest. 

It is speculated that, if a regularization is imposed to suppress the very high j contribu- 
tions, Cj 1 ' 1 can be tamed to converge to constants. More explicitly, a simple example is to 
regularize the matching condition (13 .24ft as 



n/2 

£ 

3=1/2,1, 



e -6f(j) c ( n ) 



n 



8nG-f 2 (x 



2,72 



-2k 



k sin ixc \p[ 



(3.27) 



k=l 



by introducing a monotonically increasing function f(j) and a minuscule regulating parame- 
ter 6 giving rise to a cut off for high j. With this regularization, Cj are expected to converge 
to a set of constants. 

(n) 

The peculiar feature of divergence of Ca associated with lar ge j is reminiscent of the 
infrared diver gen ce encountered in many spin-foam models (see 22|, |23| for the Ponzano- 
Regge model, [24| for the 4D BF model, and 25[] for a general review on spin-foam models). 
It is an infrared divergence, since it regards large j, namely, large lengths or areas, on the 
faces of spin foams. A renormalization procedure has been developed to divide away the 
divergence in the Ponzano-Regge model [22l . 23]; the same technique might be applied here 



to mitigate the undesired divergence (see the remark below for more comments). Another 
appealing way to get rid of the infrared divergence of the spin-foam models is to replace the 
representation theory of the group (SU(2) for the Ponzano-Regge model and SO (4) for the 
4D BF model) with that of the quantum group (SU(2) q and SO(4) q , respectively, with q 
chosen to be a root of unity). This leads to the Turaev-Viro model 2y] and the Crane- Yetter 
model 27J, respectively. On both models, it can be argued that the quantum deformation 
of the group simply corresponds to the addition of a cosmological constant in the classical 
action, the value of which is given proportionally by 1/q (see [28j for this subject). It is 
tantalizing to speculate that the quantum deformation could also be used to remove the 
divergence of , if this is indeed the case, the regulating parameter 5 can be interpreted 
as a consequence of a small but nonzero cosmological constant. 

The fact that the linear superposition of ^C'^^, if suitably regulated, can yield the 

Hamiltonian operator C'^S) bolsters the idea that it might be more natural to take into 
account all j representations when the Hamiltonian constraint is quantized, although we 
have no theory yet to derive or even suggest the values of Cj from the first principle. The 
resemblance between the divergence of Cj and the infrared divergence in spin-foam models 



6 Solving (|3.24[) . we have to deal with the linear equation Ax — b with A being anxn matrix. When n is 
large, A contains both very big and very small entries. As a consequence, the double precision used in our 
computer program loses necessary accuracy and results in noticeable round-off error. To avoid technical 
difficulty, the numerical computation is done only up to n = 15. Because of this limitation, we can only 
know the rough profile of Cj . That is, the asymptotic behavior in (|3.26p is only an outline suggested by 
the limited data of numerical computation and the literal expression should not be taken too seriously. 
Other formulae such as ~ Cj n a e n / 2 with awl are also possible. 
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FIG. 1: Coefficients cf ] scaled by (ne™/ 2 )" 1 for n = 2, 3,..., 15. {cf~ l) is a single point at 

j = 1/2, which is out of the plot.) The scaled coefficients c*- /(n e n//2 ) tend to converge to a set of 
constants c,- as n increases. 



adds one more rationale to suggest that the inclusion of higher order holonomy corrections 
in the Hamiltonian operator may come out naturally from the spin-foam models (see also 
[Til]), although the link between LQC and the spin- foam formalism is far from clear. We 
leave these issues to future research and in the second half part of this paper we will simply 
take (12.91) as our departing point to construct the corresponding quantum theory of LQC 
and study its dynamics. 



Remark: 

In the full theory of LQG, the classical Hamiltonian constraint has to be regulated before 



it is quantized. The procedure of the general classical regularization is described in [29 
Assign a partition of £ into cells □ of possibly arbitrary shapes such that every cell □ has 
linear dimension smaller than e in coordinates. For every cell □, we define edges sj and 
loops j3u (lying on the surface spanned by si and sj). Finally, fix an arbitrary chosen 
representation j of SU(2). The entire structure is denoted as R e and called a permissible 
classical regulator if 

\im C^ C \A, E) =C Euc \A,E), (3.28) 

^Ir 1 = J2 C %z (3-29) 
□ 

N(vn) 
(87rG) 2 7 3 / 2 



where 



= T^Mm E C ""^ - li) K^"h, K ,V}] (3.30) 



UK 



with fixed constants Cuk and if a similar condition also holds for the Lorentzian part. If 
we go one step further to generalize this procedure, we could remove the restriction to a 
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fixed representation of SU (2) and thus include all generic j representations. The notion of 
permissibility of the classical regulator will be slightly changed as f 1 3 . 3 j) is modified to 



cSS = j^Mm E - U %L) w Kl{ u K K , v}} (3.3i) 

1 ' 1 j,IJK 

with V'Cuk being constants. 

In the context of LQC, the naive analog to (I3.28j) would be 

^^v—iJ^V (3.32) 

for a fixed j. However, the remarkable feature of LQG that the dependence of the regulating 
parameter e disappears is no longer the case in LQC, because the diffeomorphism invariance 
is explicitly broken with homogeneity imposed. As a consequence, we should keep Ji finite 
by hand to impose the fundamental discreteness of LQG. Nevertheless, if we consider all 
generic j representations rather than choose a fixed one, the notion of permissibility can be 
restored and the analogous criterion is given by 

with 

^ = E c i n)a) ^av (3-34) 
j 

In a sense, the role of the regulating parameter e in LQG has been traded for n in LQC 
to achieve permissibility and, as a trade-off, the aforementioned divergence regarding large 
j implies that the dependence of n cannot be trivially removed. From this perspective, 
this divergence is analogous to that in the Ponzano-Regge spin-foam model, in which the 
finite triangulation (analogous to keeping ]x finite) breaks the continuous diffeomorphism 
invariance and a renormalization procedure taking into account the translation invariance 
has been developed to get rid of the divergence (23]. A regularization scheme for the case 
of LQC might also have the same rationale. 

Remark: 

The virtue of adopting (I3.1ip . instead of (13.91) . and avoiding Thiemann's trick (by scaling 
the Hamiltonian with an appropriate N) is that the resulting Hamiltonian operator for a 
given j involves only even powers of sin ftc as shown in (I3.22p . giving a systematic algorithm 
to uniquely determine the coefficients Cj via (13.241) . Had we used (13. 9p and/or taken into 
account Thiemann's trick, the resulting operator for a given j would be more involved (see 



13| for the explicit expression of ^C gr!w ) and the matching condition in parallel to (I3.24p 
would be much more complicated. This does not necessarily impede the linear superposition 
of ^Cgrav to yield Cg^l as in (12.101) . but to show the possibility is much more difficult than 
in the alternative quantization scheme we have adopted. Furthermore, if Thiemann's trick 
is used to recast the term involving the inverse triad (although not strictly necessary as 
remarked at the end of Sec. IIII Bp . correspondingly, for the matter part, the inverse volume 
operator should also be defined a la Thiemann's trick as shown in (I3.17p . which gives rise 
to both j and £ ambiguities. As well as the gravitational part, the matter part of the 
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Hamiltonian operator shall include all contributions of j. However, unlike the gravitational 
part, there is no clear sense to posit a linear combination of ^V -1 as the most natural 
choice. 



IV. LOOP QUANTUM COSMOLOGY WITH HIGHER HOLONOMY CORREC- 
TIONS 

This section presents the second half part of the paper. As the correspondence between 
higher order holonomy corrections and higher j representations of the Hamiltonian operator 
has been explored in Sec. IIII[ for the second part, we take (12.91) for granted. Its quantum 
theory of LQC will be formulated and studied by applying the techniques introduced in 
[ill ]. I n Sec. IIV A| we construct the physical Hilbert space and identify the complete set of 
Dirac observables. In Sec. IIVBI the formulation is recast in the b and x representations to 
facilitate further calculations. In Sec. IIVCI we adopt the notion of singularity resolution 
suggested in and prove that, for any arbitrary physical states and at any order n, the 
classical singularity is resolved as the matter density remains bounded from above by an 
absolute upper bound, which equals the critical density obtained in the heuristic analysis of 
[I2I ] . Finally, in Sec. IIVD} we elaborate on the bouncing scenario for the case of n = 00 (all 
orders of holonomy corrections included) and study the evolution of coherent (semiclassical) 
states. 



A. The full Hamiltonian constraint and the physical Hilbert space 

To impose the discreteness parameter p,, following the "improved" dynamics suggested 
by jl], the discreteness parameter ft is designated as 




where A 2 = A is the area gap in the full theory of LQG and A = 2^/3%^^ for the standard 
choice (but other choices are also possible) with £p\ := \J Gh being the Planck length. With 
this choice, it is convenient to introduce the new canonical variables: 



" := (Sgnp) 2^' (42a) 



IPl' (4 ' 2b) 

which, by (13.21) . satisfy the canonical relation 

{b,u} = 2h-\ (4.3) 

Also note that Jic = Xb. 

To construct the kinematical Hilbert space 7ikin of LQC, it is convenient to use the v 
representation in which states are wave functions ^f(u) and the operator measuring the 
physical volume of V is given by 

VV(v) = 27r4i7|£|*(f) = 27r4r7|f ( 4 - 4 ) 
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As before, the variable b is canonically conjugate to v. However, unlike z>, the operator b 
is no longer well-defined in LQC and should be replaced by the holonomies, which can be 
promoted to operators with the action: 



(4.5) 



The classical Hamiltonian constraint (12.91) will be taken as our departing point for the 
quantum theory. For the matter part, we adopt the ordinary Schrodinger representation in 
which is a multiplicative operator and p^ acts as a differential operator: i.e. <f>^(v,<f>) = 
4>^{y, (ft) and p^(u, <fi) = —ihd^ (v, </>). Now we are ready to promote (12.91) to the quantum 
Hamiltonian equation: 



-3^G\0\b^\v\b^(u. 



(4.6) 



where b^ is defined as 



{2k)\ 



A^2»(*!)»(2* + l) 



sin Xb 



2k+l 



(4.7) 



and a particular factor ordering has been chosen to concur with the ordering adopted in [TlJ . 
Additionally, physical states must lie in one of the irreducible representations (symmetric 
or antisymmetric) of the orientation reversal: n : v i — > n(z/) = —u, which is regarded as a 
large gauge transformation. Since there are no fermions in this model, \I/ is assumed to be 
symmetric: U.^f(u,(p) := ty(—v,(f)) = iff(y,<p). 

By taking this symmetry into account and expressing out the action of sin Xb explicitly, 
the Hamiltonian equation (14.61) is cast as 



£ 

k=0 
2n+l 



£00 U n ) 
(2k) 



2 2k (k\) 2 (2k 



(2k') 



^ 2 2k '(k'\) 2 (2k' + 1) 



k'=0 



J 1 ,J 2 =-(2n+l),-2n+l, 



(_l)(2fc+l-il)/2( 2 jfc + i)! 



x (v + hX + l 2 X) *(z/ + 2 ZiA + 2 Z 2 A) 



■i) 



(2fe'+l-i 2 )/2 



(2fc' + l)! 



(2z) 2fc,+1 ( 2fc, + 1 ~ b )! ( 2k '+i +h ) ! 



(4.8) 



which can be regarded as an "evolution equation" evolving the quantum state in the internal 

time 4>. The kinematical Hilbert space is given by 7Yki n — £ 2 (^Bohri IH~ ^ u Bohr) ® L 2 (M, d(j)). 

>00 



The factor l/|f| is prescribed in the measure to have the operator self-adjoint. For any 

order n, 0|"| is a difference operator in steps of 4A. For each e G [0, 4A), let TCj iri denote 
the subspace of the kinematical Hilbert space 7ikin with states whose support is restricted 
to the lattice points v = e + 4nA (with n G Z). Each sector 7i kin is preserved under the 
evolution equation; hence, there is superselection among these sectors under dynamics. We 
can then focus on one of the sectors 7i kin © T~&~ e if symmetry is also concerned. In this 



paper, without losing generality, we will focus on 7i kin ° for simplicity (note 7i kin ° 



oy4A 
rL kin 
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The treatment here should be straightforwardly extended for other sectors without much 
difficulty. 

The physical Hilbert space TC p hy can^be obtained by applying the group averaging pro- 



cedure to A4.8P ; the same strategies in [U| to construct the physical Hilbert space for the 
familiar n = case can be exactly followed. For generic n, again, the physical states satisfy 
the "positive frequency" square root of (14. 8ft : 



» 

and the physical scalar product is given by 



-id^{v^) = J®^(v ) <p) (4.9) 



^ 2 ) phy = ^ E i^T* 1 ^' <^)^, ^o), (4-10) 

v=An\ 

where the right side is evaluated at any fixed instant O . The physical scalar product is 
independent of the choice of 0o (i-e. conserved in internal time 0) because G)|"j (and hence 
its square root) is (positive definite and) self-adjoint with respect to (-, -) p h y given above. As 
in the case of n = 0, the state with support at v — yields infinite norm and hence does 
not belong to the physical Hilbert space. This is a consequence of quantum dynamics since 
(•, -)phy is determined by the constraint operator. 

Finally, we can define the complete set of Dirac observables. The first is the momentum 
p<j> defined as 



p^(u, 0) = -m ^ = hy/e™ 0) (4.ii) 

and the second is a 1-parameter family of the "relational" observables V\ ( f >Q defined as 

V\^{vA) = 2tt4i7 exp (iyfefy (0 - <M) WM^M, (4-12) 

which represents the volume at time 0o by first freezing the solution ^f(i/,(p) at = 0o, 
acting on it by the volume operator V defined in ( I4.4p . and evolving it from 0o to by (14.91) . 
Obviously, both p^ and V\^ preserve the solutions to (14. 9 j) and are self-adjoint with respect 
to f)4.10p since the square root of G)|™j is self-adjoint. 

To summarize, the physical Hilbert space 7i p h y (associated with the superselected sector 
H^®) consists of solutions ty(is,(f>) to ( 14. 9 j) (with support in the set {v = 4nA}) which 
are symmetric in v and have finite norm with respect to (I4.10p . A complete set of Dirac 
observables are given by p$ and V"|^ as defined in ( 14.111) and (14.12p . respectively. 



B. The b and x representations 

Since \I/(z/, 0) G TC^ has support on the "lattice" points {y = 4n\;n e Z}, the state 
in the b representation gives the wave function 

:= (W = E (WW*) = E ^ 6 %,0) (4-13) 

v=An\ v=An\ 
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as the Fourier transform of ^>{y,<f). It follows that ^{b + 7r/A, 0) = \I/(&, 0) is a periodic 
function and the inverse Fourier transform reads as 

f := (i/ltf) = - f 1 db(u\b)(bm = - f ' dbe-$ ub ty(b,(j)). (4.14) 

* Jo ' ' 7T Jo 

Additionally, the symmetry requirement \l/(z/, 0) = \1>(— z/, 0) is equivalent to ^(6,0) = 
6, 0). In 6 representation, z> and sinA6 act as 

£*(M) = -2ij^(M) (4-15) 

and 

slmA6*(6,0) = (sinA6)^(6,0). (4.16) 

Following the treatment used in set 0) := A/(7rz/)\l>(z/, 0) and let xip,(p) := 
^„ = 4 n ^ e ivh ' 2 x{y, 0) be the Fourier transform of x(v,(j)). Then, on x(6, 0), the Hamiltonian 
constraint (14. 8 p becomes 

dlxiPA) = 12vrG (&i n) <9 fc ) 2 X (&,0) =: -e$x(M). (4.17) 

Also note that the symmetry requirement \&(— z/, 0) = \l>(z/, 0) is now translated to 

X(~v, 4>) = -Xiy, 4>) and 0) = ~x0>, 0). 

In the b representation, the physical scalar product (I4.10p takes the form 

ptt/X 

(Xl>X2)phy=/ ^Xl(&,0o)|i > |X2(&,0o) ) (4.18) 

where \v\ is the positive square root of the self-adjoint operator u 2 = (—2idb) 2 on Wj^jf = 
L 2 ([0, 7r/A), db). 7 The Dirac observables defined in (14. lip and (14.121) now read as 

M b ></>) _ fc faJn) 
dcj) 



PM b i <t>) = -^- 1 kr z = V e X ^) ( 4 - 19 ) 



7 To arrive at (|4.18[) from (|4. 10|) . we have to apply the formula for the Fourier transform of a Dirac comb: 



e -i2n7r/T 



n— — oo 



which leads to 

oo 



A ^ V A 

t/— 4nA k— — oo 

= ~5(b-b'), if 6,6' e [0,tt/A). 
A 



This equation can also be used to show 



*(&, 0) = - (-2i^) X (6, 0) = y »x(b, <f>). 
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and 

V% oX (&,0) = 27r4 l7 expfzy^^-^,)) |i)|x(Mo)- (4-20) 
To simplify the Hamiltonian constraint further, let us define the new variable x 

f b db' 

x= - — for b G [0,tt/A), (4.21) 

such that \/Y1txG b^'dh = d x and x G (—00,00). For generic 6, the corresponding new 
variable x is then defined as x(b) = x(bo) with 60 — b — nir/\ G [0, 7r/A) for some n G Z. 8 
Particularly, for the case of n — 0, we have = sin(A6)/A and 

I= 71b ln(ta,1 T)- < 4 ' 22 ' 

For the case of n = 00, it follows from (14.71) that bh : = b h °°^ is a triangle wave with period 
27r/A and yields 

La u > 2A — u — a ' 
when restricted to [0,7r/A]; 9 correspondingly, the new variable x reads as 

2A 



In ^b, for < b < 



x _ ^ mo ^2A /tt _ ft \ f jr. < 6 < 2 * ( 4 -24) 

Tl2^G m 7T VA ^ ' iUI 2A — — A ' 

For other cases of n, the expression of x as a function of b is much more complicated. 
In terms of the new variable x, (14.171) becomes the familiar Klein-Gordon equation: 

dlx(x, 0) = d 2 xX (x, 0) =: -G)g X (x, 0)- (4-25) 

Again, the physical states can be taken to be positive frequency solutions to (14.251) . i.e., 
solutions to 

-ifyx(z,0) = ij$$x(x,<f>). (4.26) 
If the initial data at the instant <fi = 0o are given by 

x (x,0 o ) = -=/ dke~ ik *x(k), (4-27) 



2,71 



8 By this periodic definition, the 7r/A-periodic function (f>) is well posed as x{. x i4>) '■— x(b,<fi) with the 
new variable x € (— 00,00). For generic b, we have \J\2-kG \b h \db = d x and the symmetry requirement 
x(— b, (f>) — — x{b, 4>) now reads as x(— x, <f>) — — x( x i <t>) m terms of the new variable x. 

9 Note that b^' (bh included) is 27r/A-periodic and antisymmetric under b — >■ —b, but nevertheless the 
operator {b^ db) 2 and its positive square root are 7r/A-periodic and symmetric under b — ► —b. Thus, 
when acting on x(b,4>)> the operator (jj^db) 2 and its positive square root do not alter the 7r/A-period 
and antisymmetry of %(6, <p). In particular, bh is a triangle wave with period 2tt/X and agrees with b only 
for b e [— 7r/(2A),7r/(2A)]. The distinction between bh and b is an essential consequence of the intrinsic 
discreteness of LQC. That is, for any € 7ikin, ^(fi, <fi) is an almost periodic function of b, which becomes 
periodic if restricted to HL . 
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then the physical state as the solution of (I4.26P takes the form 



oo 




I f I f 

— I dke- ik ^ +x) e lk4,0 x(k) + -= / dke ik ^- x) e- lk *°x(k) 
'2n J -oo v27r jn 



= XL(<l> + x) + XR((t>-x)- (4.28) 

As usual, the solutions can be decomposed into the left and right moving modes. The group 
averaging procedure implies that the physical product is the standard scalar product for the 
Klein-Gordon equation: 



(Xi,X2) P hy = -i dx[xi(x,<P)dtX2(x,<f>)-(dtXi(x,<p))x2(x,<i>)], 

— oo 
oo 

dk\k\xx{k)x2{k), (4.29) 



which is conserved in the internal time 0. Using d^XL(<P+ x ) — 9 x xl(4 ) + x ) an d d^XR^—x) — 
—d x xii{4>—x), it can be shown that the left and right moving sectors are mutually orthogonal, 
i.e., (xiL,X2i?)phy = 0. This follows 



■DO 



(Xi,X2) P hy = -2i / dx [xil((P + x) d x X2L((j) + x) - Xm(4> - x ) d x X2n(4> - x)} . (4.30) 

J — oo 

Noting that the operator v = —2idb = —2i{12nG)~ l l 2 b^ d x is positive definite on the left 
moving sector and negative definite on the right moving sector, we then have 



(Xi,X2) P hy = / dxxi(x,(f))\-2id x \x2(x, 

-oo 

i-tt/X _____ 



db xi (x(b) , 0) | -2id b | X 2 (x(b) , 0) , (4.31) 

where \—2id x \ and \—2idb\ denote the positive parts of the self-adjoint operators —2id x and 
—2idh on L 2 (R, dy) and L 2 ([0, tt/X), db), respectively. This equation agrees with (I4.18p . 
By ( I4.19P and (14.201) . in the x representation, the Dirac observables act as 

p^xixA) = -iftd^xix,^) = -ihd x XL(4> + x) + ihd x XR(4> - x) (4.32) 

and 

VUoX(x, 0) = 2*4,7 eV 5 ^ \v\ x (x, O ) (4.33) 
= 2*4,7 eV 5 !^*' [v X l(<Po + x)- v X r((Po ~ x)] 



2i 2?r£ P' 7 e ^^-M 
2*4,7 



b^\x) 1 <9 2; xl(0o + x)- b£Xx) l d x XR(<Po - x) 



-2% 



6j t n) (a; + 0-0o) O x xl{<P + x) -b ( ^\x-4>+(f)o) d x XR{<f) - x) 
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Taking into account the symmetry requirement \{~ x -> 



-x(x, 0), we have 



X (x, 0) = + x) + -x) = -j= + x) - F{<f> - x)) , (4.34) 

i.e. Xl(x) = —Xr(x) = F(x)/V2. Equivalently, this is to say x(— k) = —x{k) in (14.281) . 
The full information in any physical state xi x i 0) is contained in F, and we can conveniently 
describe 7i p h y in terms of the positive frequency, left moving solutions F(<f> + x) (or, equiva- 
lently, the right moving solutions F((p — x)), which are free of any symmetry requirement. 10 
In terms of F, the physical scalar product (I4.30P is simply written as 



(Xl,X2, 



phy 



-2i 



dxFAx) d x F 2 (x). 



(4.35) 



Equations ( 14.32(1 and ( 14.33(1 then give the matrix elements of the Dirac operators as 



phy 



-2H 



dx [xil(<P + x) dlx2L((j) + x) + xir(<P ~ x) d 2 x X2R{4> - x)] 

) 

/CO 
dx[d x x 1L ((f) + x) d x X2L(4> + x) + d x xm{4> ~ x) d x X2n(<f) - x)] 
-00 

POO 

2h / dx d x F(x) d x F 2 (x) (4.36) 



and 



<t>oX2) phy 



4(2tt4 i7 ) 
VV2^G J 

2(2tt4 i7 ) 



VUttG J 

Similarly, for later use, we can also derive 



dx d x XiL{<Po + x) b^\x) d x X2L{<Po + x) 

+ dxXmi^o ~ x) b^\x) d^R^o - x) 
dxd x F\(x 



(4.37) 



b^\x - O ) + btXh - x)' 1 } d x F 2 (x) 



{xupI X2) 



phy 



.Xi,^X2) P hy = -2ift 2 / dxd x F 1 (x) d 2 x F 2 (x 



(4.38) 



and 



(XU (V r Uo) 2 X2)phy = (VlfoXl, V r UoX2)phy 

-8mGh 2 ~f 2 r °° 



dx ( b^\x 



or%F(x) ) d x (bPix - <p )- l d x F(x) 



(4.39) 



Note that although F{x) is free of any symmetry requirement, F(x) is not an arbitrary function but the 
function such that F((f> + x) is the left moving solution (or, equivalently, F(cj> — x) is the right moving 
solution). The general solution is given by (|4.56|) . 
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C. Singularity resolution and the absolute upper bound on matter density 



We have noted for (I4.10p that the state with support at v — does not belong to the 
physical Hilbert space 7i p hy This by itself, however, does not mean that the singularity is 



111 ] for more comments on weaker notions of singularity resolution.' 



avoided. (See 

In fact, it is nontrivial to obtain necessary and sufficient conditions to characterize the 
occurrence of a singularity even in classical general relativity, and thus a generally applicable 
and satisfactory notion of "singularity resolution" is not available in quantum gravity. In 
simple situations, however, the notion of a singularity is unambiguous in the classical theory 
and a satisfactory notion of singularity resolution is possible in the quantum theory. A 
sensible notion is proposed in [ill ]: The singularity is said to be resolved in the quantum 
theory if the physical Hilbert space and a complete family of Dirac observables are available 
and the expectation values of these observables remain finite in the regime in which they 
become classically singular. In the model studied in this paper, the matter density diverges 
at the classical singularity. The singularity is then considered as resolved in the quantum 
theory if we can show that the expectation values of a complete set of Dirac observables 
including the matter density remain finite. 11 

For the case of n — 0, it has been shown in 



111 ] that the matter density is bounded by 



an upper bound as far as two natural notions representing the mean value of the matter 
density are considered. We will show that the same is also true for the cases of generic n, 
while the upper bound is modified by a numerical factor involving n. 

First, let us define the mean value of the matter density (as a function of the internal 
time 6) as 



{x,h x) 



phy 



(4.40) 

2(x,VUx) p 2 hy 

and use it to represent the measurement of matter density. From (14.361) and (14.371) , it follows 



JZdx\d x F(x) 



Noting that (T4T7]) gives 



J_" dx \d x F(x] 



i 



i\i 2 



X 



(4.41) 



(2k)\ 



Xf^2 2k (k\) 2 (2k + 1) 



we then have 



M0) < 



f°° dx\d x F\ 

J— OO 1 x 1 



2*G>f (2\r n 1 ) 2 r oo dx\d x F\ 2 
where the Planckian density is defined in ( 12.51) . 



3n 

T' 



3 SnPPl 



(4.42) 



p (n) 
rsup' 



(4.43) 



Note that the exclusion of the state with support at v — does not imply that the expectation values of 
the matter density would be bounded from above on the physical Hilbert space H p h y , because does 
not have an upper bound on 7i p hy 
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This shows that the mean density p</,(0) is bounded from above by an absolute upper 
bound piup, which exactly coincides with the critical density p^} t obtained at the level of 
heuristic effective dynamics Also note that 



-JPPI — Psup ^ Tsup ^ Ps 



37T 2 



sup 



-ppl- 



(4.44) 



(n) 

That is, the upper bound ps Up increases with n, but remains finite in the Planckian regime 
even if all orders of holonomy corrections are included. 

Second, we can use the expectation value of the density operator in place of p^ used 
above to represent the mean density. To have a self-adjoint operator, let us define the 
density operator at the instant as 



(4.45) 



Analogous to (14. 33ft . the inverse volume operator V\^ at the instant 0o is defined as 



(VU-'xix^) := (2tt4 i7 ) 



-1 



(x) 



*°) U,-l 



xO,0o) 



(4.46) 



^ (bfxL^o + x))- d- 1 ( C j X^(0o - x) 



where is the inverse of |-0 1 , 12 and d x 1 is the inverse of d x . 13 Consequently, by (14.321) . it 
follows 



4vr4 l7 



hi hX x )XL{<Pv + x)+ &J™(ar)Xfl(0o - x) 



(n)/ 



47r4i7 

By (14.351) . we then have 



b ( ^\x+(f)-(f) )xL(4> + x) + b^\x-(j)+(f) )xR(4> ~ x) 



(4.47) 



phy 



-62 



47rG 7 2 

— 6? 
47rG 7 2 



-oo 
oo 



dx ( ^(X + - 0o ) ^(0 + X) ) d x ( fiJ^X + - 00 ) ^(0 + X) 



( b ( Z\x - O ) F(x) ) d x ( 6™^ - 0o) F(x) ) . 



(4.48) 



12 Note that |i> _1 | is not densely defined on 7ikin, but nevertheless it is well defined on H p hy since the state 
with support at v = is excluded on 7i p hy 

13 More precisely, when acting on x(x,4>), the operator d~ l is given by 

d- l x {x, $ = -= \ dk — e-^+Am-M m . 
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Therefore, the expectation value of the density operator yields 



,. fc(i > U)" 1 ^(v'U)- 1 x)phy 

M4,)) : = 2( X , x)pby ' 

3 (b^Xx - (p)x, b ( h\x - 0)x) P hy 
8vrG7 2 (x, x) P hy 

<T 3 (A 1 $ n ) 2 (X> X)phy _ Q ^2 n _. An) (a aq.\ 

- a n 2 1 \ - 3^ n p P i -. p y su > (4.49) 

87TG7 2 (X, X)phy 

where (I4.42p has been used again. This shows that, as the second notion of the mean density 
is regarded, the matter density is again bounded from above by the same upper bound psup- 
To summarize, both p^ and (fy) can be used to represent the physical measurement of 
the matter density, and both of them are bounded from above by the same absolute upper 
bound piup, which equals the critical density p^t obtained at the level of heuristic effective 
dynamics. Therefore, we have shown that, at any order n, the singularity is resolved as 
far as the above two notions of the matter density measurement are concerned. Moreover, 
the singularity resolution holds for any arbitrary physical states, not only restricted to the 
states which are semiclassical at late times. 



D. Quantum bounce and the bouncing scenario 

We have shown that the matter density is bounded from above in both senses of the 
matter density measurement. For the case of n = 0, it can further be shown that, for any 
arbitrary (not necessarily semiclassical) states, the expectation value of the volume operator 
(V\^) follows the bouncing trajectory and the matter density at the bouncing epoch can 
come arbitrarily close to piup ^ In the following, we will prove that the same conclusion 



can also be drawn for the case of n = oo. (For the case of generic n, we expect the same 
but the algebra is much more complicated to give a rigorous proof.) 
Equations (14331 and KMf yield 



b h (x) := bt\x) = \ fOT — < - < 0' (4.50) 

^e-^G*, for 0<x<oo. V ; 



Taking this into (I4.37P then gives 

(VU) = (x,v\<px) P hy 



where 



/(p poo 
dx\d x F(x)\ 2 e- VT ^ {x -^ + / dxldxFix^e^^^ 
, 00 J(j> 

V+((/))e vmu<l ' + V^(((>)e- y/mu ' 1 ', (4.51) 



V ty) = r dx\d x F(x)\ 2 e-^ x =: t dxv + (x), (4.52a) 

= 167^ r dx \ d F r x)l 2 e+ VmG,^.^3^ HdxvJx). (4.52b) 
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Similarly, taking (14.501) into (I4.39P gives 



((^U) 2 ) = (x,(^U) 2 x)p hy = ^ + (0) 



<t> + W_((p) e 



-2Vl2nG<, 



(4.53) 



where 



32 GfrVA 2 



W + (0) 



dx- 



1 



3vr j_ 00 i 

2 .,U 2 



d x 



-sJ\2-kGx 



d x F(x) 



W_(0) 



32 (4i7 A) 2 

3ttG 
32G^VA 



— oo 

2_ 2 \ 2 /-oo ^ 

, dx— 

37T ./^ Z 



(4.54a) 



<9,: 



32 (4 l7 A) s 
3vrG 



dxit;_(x). 



(4.54b) 



As long as F(x) is smooth, the trajectories of (V^), ((^|^) 2 ), and thus (AV|«^) 2 : = 
((Vltj)) 2 ) — (V\,f,) are all smooth functions of 0, despite the kink of the function bh(x) at x = 0. 
This tells us that the abrupt kink of the solution obtained in 12| at the level of heuristic 
effective dynamics is only an artifact, which is smeared by the quantum fluctuations at the 
level of quantum theory. Furthermore, for the physical states which are highly semiclassical, 
the integrand in (I4.52p and (the real part of ) the integrand in ( 14. 54[) are appreciable only 
for some localized intervals and diminish rapidly away from the intervals; i.e. v±(x) ~ 
if x ^ (x~,x+) and (the real part of) w±(x) ~ if x ^ (x~,x+) for some x^ and x^. 
Consequently, we have V + (<p > x+) ~ V+(oo) =: V + and V_(0 > x+) ~ V_(oo) = 0; 
r_(0 < x-) ~ VL(-oo) =: V- and V + (<p < x~) ~ 0; W + (<f> > x+) ~ W + {oo) =: W + and 
W_((f) > x+) ~ 0; as well as W_(0 < x~) ~ W_(-oo) =: W_ and W+(0 < x~) ~ 0. This 
leads to 



(VU) 



y Vl2^G^ for ^ > x j 



-V12ttG<; 



for 



for 
for 



< x; 



(4.55a) 
(4.55b) 



with V± and W± being constants. 14 Therefore, in the distant future and past (i.e. > 
max(x+,x+) or < min(x~,x~), respectively), (V^) follows the classical trajectory with 
constant relative uncertainty spread (i.e. (AV\^) 2 / (Vl^) 2 ^ (W±—V 2 )/V 2 ). In other words, 
(I4.55P gives exactly the same behavior as the WDW theory in the far future and past. Thus, 
the physical solution (I4.5ip gives rise to the bouncing scenario in which two WDW solutions 
(expanding and contracting) are bridged together through a transition phase of the quantum 
bounce. Additionally, the more the integrands v±(x) and w± e (x) are localized, the more 
abrupt the transition phase is. In principle, the transition phase can be arbitrarily brief if 
v±(x) and w± e (x) are extremely sharp. 



14 Note that while v±(x) is real, w±(x) is complex in general. However, the contributions of the imaginary 
part of w±(x) should exactly cancel out in (|4.53|) since ((V^) 2 ) is real. 
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To understand the dynamics in more detail, let us study the evolution of Dirac observables 
for coherent (semiclassical) states explicitly. According to (I4.28[) and (I4.34p . the general 
solution of F(x) is given by 



F(x) = V2 X l(x) = -V2 X r(x) 



/ 27T 



dke~ lkx F{k) 



dke +ikx F(k) 



(4.56) 



where, without losing generality, the offset constant <f) in (I4.28[) is set to <f) = 0. 15 To 
construct a coherent state, we choose the Fourier amplitudes to be 



(2.7)2 



F(k) = -V2x(k) = ^ZXI 



e (2^p 



V / 2(27ro- 2 ) 1 / 4 VfcT 



for k > 0, 



for fc < 0, 



(4.57) 



where F(k) = —F(—k) and -F(fc) takes the form of a Gaussian distribution centered at k* 
(k* > 0) with a representing the width. In case of highly semiclassical states where the 
Gaussian distribution is very sharp with a <C k*, we can accurately approximate (I4.56P to 



F(x) 



dk e 



-ikx 



a 



-cr 2 x 2 +ik t x 



'2-irJ-oc v / 2(2tt ( t 2 ) 1 /4 v ^; (2tt) 1 /4^ 

By (14.351) . (I4.36P and (I4.38p . it is straightforward to show that (x,x)ph y ~ 1 and 

(P4>) = (x,P<t>x) P hy ~ UK (l + (cr/K) 2 ) , 
(#> s (x,pJx) P h y ^^ 2 (l + 3( ( T/^) 2 ), 
(A^) 2 = (p}> ~ W « ^ 2 (1 - (a/fc,) 2 ) . 



(4.58) 



(4.59a) 
(4.59b) 
(4.59c) 



That is, the coherent state given by (14.581) is normalized and the parameters k* and a are 
associated with the expectation value and the uncertainty spread of p^. With a <C fc*, ( I4.59P 
then yields (p^,) m hk* and (Ap^) 2 rs h 2 a 2 . 

Given with (I4.58p . the integrand in (14.521) reads as 



v ± (x) := \d x F(x)\ 2 e TVT2 ^> 



a 



2nk* 



(4crV + k 2 ) e -^ 2 - 2 TVT^Gx 



(4.60) 



and the integrand in ( 14.54ft reads as 



w±{x) 



eT Vl^Gx Q xF ^ 



1_ g-20- 2 x 2 ^-2\/\2-kG x 



2i\k* 



d x e^" M ^ x d x F{x 
K {k 2 + 2a 2 + 4aV) 



+ i 



--: wf(x)+iw l ™(x). 



kl [2a 2 x± V12tcG) + 4 a 4 x -1 ± VUnG x + 2 a 2 x 2 



(4.61) 
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As we will see shortly, this choice will set the epoch of the quantum bounce at = 0. 
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Note that for the imaginary part of w±(x) yields 



dxw l ™(x) 



oo 
oo 



dxw^(x) 



2y2nk^ 
a 



+ 



2V2ttK 



(4.62a) 
(4.62b) 



The contributions of the imaginary part of w±(x) cancel out exactly in (14.53p as expected. 
We can simply ignore w^ l (x) and replace w±(x) with the real part w± e (x). 

Equations (I4.60P and (14.6 ip show that both v±(x) and w± e (x) are localized bumps with 
widths which decrease as a increases. This confirms what we have just discussed prior to 
(I4.55p . For a given k*, the bigger a is, the narrower the bumps are and thus the briefer the 
transition phase of the quantum bounce is, as can been seen in Fig. [21 

Substituting (14.601) into ( 14.521) and (the real part of) (14.61j) into (14.541) . we have 



V±(<j>) 



V2 



a 



+ e ^ {kl + a 2 + 3ttG) 



1 + erf 



V3^G ± 2a 2 c, 



V2 



a 



(4.63) 



w ± {<p) 



(4i7A)^ 
3vrG 



4a 
n/2tt 



+ (^ 2 + 3a 2 + 127rG) 



1 + erf 



V6^G± V2 cr c 



a 



, (4.64) 



which give the evolutions of (V|^) and (AV|^) by (14.511) and (14. 53ft . The evolutions of (V|^) 
and (AV U) 2 /(VU) 2 are depicted in Fig. El 16 

Particularly, we are interested in the limiting constants: 



V+ := V±(±oo 



4 ( 7 4 L A) e 



3ttG 
2a' 2 



W± := W±(±oo) 



y/'SnGk* 
16( 7 4 1 A) 2 e^ 



3ttG 



(&* + cr 2 + 3ttG) , 
(A; 2 + 3a 2 + 12ttG) 



(4.65) 



(4.66) 



which are used to compute the constant relative uncertainty spread in the WD W regime in 
terms of k* and a: 



lim 

i-±oo (y\)2 



(w ± -v£)/v z 



e^k 2 (k 2 + 3a 2 + 127rG) 
{kl + a 2 + 3vrG) 2 



37TG 37rG r — 

~ — 1 « — — , if > cx > V7rG. 



2 ' 



(4.67a) 
(4.67b) 



To exaggerate the quantum effect, in the figure, we choose fc* 3> v irG and fc* 3> cr but not cr ^> VttG. 
As a caveat, the limiting constant of the relative uncertainty spread is given by (TQ7a|) instead of (|4.67b|) 
and it does not necessarily decreases as a increases as suggested by (|4.67b|) . 
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Finally, to know the epoch of the quantum bounce, let us compute 



d^VU) = V^)e^t> + VL{(l>) 



Wl2nG(, 



* + VT^G (V+(0) e VT ^^ - V-(<f>) e 



Vl2nG< 



+Vl2^G (V+(0) ef^* - V^cj)) e ~ y/m3 ^ 



(4.68) 



It is easy to show that d^(V\^) = if and only if = 0. Therefore, <\> = is the epoch of the 
quantum bounce and the expectation value of the volume at the quantum bounce is given 
by 



bounce 



(VVo) = V+( 



PI 



Act 



SttG 

e 2o^~ 



2irk* 



O) + y_(0 = O) 

(kl + a 2 + 37rG) 



1 + 



6G & 

(k 2 + a 2 ) } if K, a > VttG 



1 + erf 



V2a 



VUttGK 

By (I4.59ap . the matter density at the bouncing epoch is given by 

(p^) 2 3n 2 e~^ 



(4.69a) 
(4.69b) 



(do) 
/'bounce 



2V 2 

bounce 



3ttG 



Ppi 



(4.70a) 
(4.70b) 



where psup is given in ( 14.44p . Therefore, for semiclassical states (K ^> a ^> \/irG), p^ 
is smaller than but very close to pi^J . The consideration with coherent states shows that 
on H phy , p^ nce can c 
case of n = in [ll| . 



(CO) 



on Hphy, Pbounce can come arbitrarily close to pi^p . A similar result has been shown for the 



Remark: 

One of the virtues to adopt the improved scheme as in (14.1 j) is that the resulting dynamics 
is independent of the elementary cell V at the level of heuristic effective dynamics, hence 
giving the correct semiclassical behavior (see lOj for more details). In quantum theory of 
LQC, however, the invariance under the different choice of V is no longer exact but respected 
only in the semiclassical regime. This can be seen from (14.631) and (14.641) . where the constant 
\JtxG arises in the company of and/or a inside the parentheses. Since both k* and a scale 
linearly with respect to the size of V while V ttG is simply constant, the appearance of 
V txG breaks down the independence of V, which now holds only if k*,a ^> V ttG. As a 
consequence, the matter density at the bouncing epoch is slightly dependent on a (and thus 
on V as well) according to (14.701) . It is a common phenomenon that a quantum system reacts 
to macroscopic scales introduced by boundary conditions; the breakdown of the scaling 
invariance is reminiscent of the well-known "conformal anomaly" as a "soft" breaking of 
conformal symmetry. Also note that the breaking of independence of V is generic for any 
order n of holonomy corrections. 
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FIG. 2: The evolution of coherent states with fc* = 2. x 10 3 v / 8vrG and (i) a = O.Olfc* , (ii) a = 0.02/c* , 
and (iii) a = 0.03/c* . (a) The expectation value of volume (V\j,) as a function of (f> in the logarithmic 
scale. Two classical trajectories (straight lines in the logarithmic scale) are bridged through the 
quantum bounce. The value of Vbounce decreases as a increases as suggested in (|4.69p . The value 
that corresponds to pl^p is indicated by the dashed line and sets the absolute lower bound for 
^bounce as suggested in (14.70p . (b) Relative uncertainty spread (AV\,p) 2 / (Vl^) 2 , which asymptotes 
to the constant given by (|4.67ap . Notice that as a is larger, the transition phase (during which the 
relative uncertainty spread deviates from the asymptotic constant) is briefer, in agreement with 
the discussion in the paragraph after (|4.62p . 



Remark: 

Recall that in the v representation, the Hamiltonian equation (14 .8p with higher order 
holonomy corrections gives rise to a higher order difference equation. It has been argued in 
(30l | that, if the higher order difference equation admits solutions with growing amplitudes, 
the difference equation is not locally stable. Furthermore, the analysis of [15] also sug- 
gests the existence of spurious solutions in LQG for the Hamiltonian constraint in j > 1/2 
representations. It was shown in [l3] for LQC that the Hamiltonian constraint operator 
in j > 1/2 representations indeed allows spurious solutions which eventually break down 
the semiclassicality in the large scale. The ill-behaved spurious solutions suggested by the 
earlier studies call into question the validity of the higher j quantization. However, our 
investigation in this subsection shows that (at least for the case of n = oo) the expectation 
values of Dirac observables are well behaved and follow the WDW trajectories in the large 
scale, given that F(x) is a physical state. In other words, if the physical solution comes 
close to the WDW solution, it will continue following the WDW trajectory in the larger 
scale. Apparently, the problem of spurious solutions is gone in our analysis; two different 
explanations might explain the avoidance of spurious states. The first explanation is that 
the spurious solutions come out only in the quantization of a specific j but are suppressed 
if all j representations are included to match the desirable form of (I2.9p . which is motivated 
to yield better semiclassical behavior. The second possibility is that spurious solutions do 
exist in the kinematical Hilbert space but have vanishing or infinite physical norms in the 
physical Hilbert space 7i p h y and thus are excluded in 7i p h y . It is not clear whether the 
quantum evolution is still free of spurious solutions if the higher order corrections on the 
inverse volume operator as in (I3.17P are also taken into account. More studies are awaited 
for the issues of local stability and spurious solutions. 
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V. SUMMARY AND DISCUSSION 



By applying the techniques introduced in [111 ] . we have successfully formulated the quan- 
tum theory of LQC for the k = FRW model with the extension of higher order holonomy 
corrections. 

At any arbitrary order n of holonomy corrections, the physical Hilbert space is rigorously 
constructed and a complete family of Dirac observables are identified. Two natural notions 
representing the measurement of the matter density, as defined in (I4.40p and (14 .49 p . are 
used to compute the expectation value of the matter density. As far as both notions are 
concerned, for any arbitrary physical states (not only restricted to the states which are 
semiclassical at late times), the matter density remains bounded from above by an absolute 
upper bound piup as shown in (14.431) and (I4.49j) . Therefore, following the same notion of 
singularity resolution suggested in the classical singularity is said to be resolved in the 
quantum theory. We thus extend the key results of to the case of generic n. Furthermore, 
the upper bound piup given by (I4.43P increases with n but remains finite in the Planckian 
regime even for n = oo. The value of piup are exactly the same as that of the critical density 
p£rit obtained in the heuristic analysis of 

Particularly, for the case of n = oo, it is proved that the expectation value of the volume 
operator (V"|</,) gives rise to the bouncing scenario in which two WDW solutions are bridged 
together through a transition phase of the quantum bounce. Given that F(x) is smooth, the 
trajectories of (V^), ((V|</>) 2 ), and so on are all smooth functions of 0, despite the kink of 
the function bh(x) = a£° (x) at x = 0. Furthermore, the detailed analysis for the coherent 
states shows that the matter density at the bouncing epoch Pbounce ^ s sma ller than but can 
come arbitrarily close to psup . 

On the other hand, we also explore the idea that the higher order holonomy correc- 
tions can be interpreted as generic j representations for holonomies in the Hamiltonian 
constraint operator. We demonstrate that it might be possible to have a linear superpo- 
sition of the Hamiltonian operators in generic j representations match the well-motivated 
operator Cg r av°°' ) • However, the coefficients in (I3.24p diverge as n — > oo and need to 
be regulated. The peculiar feature of the divergence of cf 3 is reminiscent of the infrared 



divergence encountered in many spin-foam models [22|, [23|, [24J and the regularization may 
correspond to a nonzero cosmological constant [28[ . We hope that our study inspires further 
research on the issues of j ambiguity by exploring the link between LQC and the spin-foam 
formalism. For example, the problem of finding a crossing- symmetric linear combination of 
the Hamiltonian operators in LQG as proposed in |l4j could be related to that of finding 
the linear sum of Cgral to match the operator Cg^ 00 ^ in LQC. 

Finally, it should be remarked that the problem of spurious states associated with the 
higher j representation as suggested in 13_|, ll5| does not seem to happen (at least for the 
case of n = oo) in our model. However, we do not know whether it is because the spurious 
solutions are suppressed in the kinematical Hilbert space when all j representations are 
properly summed or because they have zero or infinite physical norms and thus are excluded 
in the physical Hilbert space. 
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